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ABSTRACT 

We present the discovery of nonradial pulsations in a hot subdwarf B star based on 30.5 days 
of nearly continuous time-series photometry using the Kepler spacecraft. KIC 010139564 is 
found to be a short-period pulsator of the V361 Hya (EC 14026) class with more than 10 
independent pulsation modes whose periods range from 130 to 190 seconds. It also shows 
one periodicity at a period of 3165 seconds. If this periodicity is a high order g-mode, then 
this star may be the hottest member of the hybrid DWLyn stars. In addition to the resolved 
pulsation frequencies, additional periodic variations in the light curve suggest that a significant 
number of additional pulsation frequencies may be present. The long duration of the run, the 
extremely high duty cycle, and the well-behaved noise properties allow us to explore the 
stability of the periodic variations, and to place strong constraints on how many of them are 
independent stellar oscillation modes. We find that most of the identified periodicities are 
indeed stable in phase and amplitude, suggesting a rotation period of 2-3 weeks for this star, 
but further observations are needed to confirm this suspicion. 

Key words: Stars: oscillations - stars: variables - Stars: subdwarfs 



1 INTRODUCTION 

Subdwarf B (sdB) stars are highly evolved low mass stars that have 
survived core helium ignition (presumably through the core helium 
flash) and now popu late the extreme blue end of the horizontal 
branch. Heber (2009) reviews the general properties of these stars. 
Generally, they are expected to have a m ass of ~0-5 Mq. Their hot 
surfac es (T eff « 22,000K to 40,000K - lHebeJll984 Saffe r et al.l 
1 19941) indi cate that they retain a thin (< 10 _3 M Q ) hydrogen sur- 
face layer dDorman et aO 19931) . While the general trends of mass 



and T c fl on the horizontal branch are well understood, the origin of 
the hot sdB stars remains unsettled. One plausible mechanism by 
which stars lose all but a small fraction of their hydroge n envelope 
has been att ributed to binary evolution and mass transfer faan et all 
l2002Ll2003h . 

The discovery of multiperiodic pulsating sdB stars 
( Kilk enny etal] 1 19971) opened up these stars to asteroseismic 
probing, promising to give us a clearer picture of their current 
internal structure and, thereby, their origins. These pulsators, now 
known as V361 Hya stars (and as EC 14026 stars, but often referred 
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to as short period sdBV stars), show multiperiodic pulsations with 
periods ranging from 1.5 to 10 minutes, with most in the shorter 
period part of that range. They typically have pulsation amplitudes 
(A I / 1) of 1 percent and le ss, with detailed st udies revealing a few 
to dozens of frequencies dReedetal.H2007al) . These periodicities 

represent nonradial p-modes of low order. 

So on aft er their discovery, theoretical work bv lCharpinet et al.l 

J1996L l200lh and collaborators identified the pulsation driv- 
ing mechanism as cyclical ionization of iron. Furthermore, they 
showed that for driving to occur, radiative levitation needed to en- 
hance the abundance of iron within the driving region. Their models 
demonstrate that equilibrium computation of diffusion effects can 
indeed produce the iron enhancement needed for driving stars at the 
effective tempe r atures (and gravities) of the V36 1 Hya stars; see 
Charpinet et al. (2001) and Fontaine et al. (2006) for an overview 
of driving in sdB stars. 

Th ere is a second, rela ted class of pulsating sdB stars, discov- 
ered bv lGreen et alj |2003). These stars, known as the V1093Her 
stars (or PG 1716 stars, or long period sdBV stars), show longer pe- 
riod pulsations ranging from 45 to 120 minutes or more, and typi- 
cally have semi-amplitudes lower than 0.2 per cent dOstensen et al.l 
2009). They are cooler than the V361 Hya stars, though there 
is some overlap, an d they are most likely g-mode pulsators 
dFontaine et alj|2003h . A third type are hybrid (DW Ly n) pulsators 
which show both long and short period variations ( Bar an et"al] 
2005; Sc huhetalj|2006l) . In most cases, the DWLyn stars show 
much larger pulsation amplitude in the higher-frequency modes, 
by factors ranging from 5 to 20 or more. 

As sdBV stars often pulsate with tens of frequencies, with 
each frequency probing a slightly different region of the interior, 
successful application of asteroseismology has great potential for 
probing the interiors of these stars. Asteroseismology of pulsat- 
ing sdB stars can provide determinations of their global properties 
such as mass and radius. Matching the pulsation frequencies with 
model values can further constrain their internal structure in terms 
of the thickness of their surface hydrogen layer and, perhaps, the 
size of the heliu m-burning core dCharpinet et al.l2007l) and internal 
rotation profiles dKawaler & Hostlerl2005l ; IVan Grootel et al.l2008l ; 
ICharpinet et"ail20 08). The mere presence of pulsations allows test- 
ing models of radiative levitation and diffusion in these stars as well 
dCharpinet etal] [2009). Taken together with spectroscopic analy- 
sis, the pulsating sdBV stars may indeed provide the long-sought 
window into these stars that can reveal thei r origins as either sin- 
gle stars, or the produ ct of binary evolution dHan et all2002ll2003l ; 
iHu et al j2003,l2009l) . 

Continuing challenges for asteroseismology lie in the stability 
of the observed periodicities, and in the identification of the degree 
and order of the observed frequencies. They are generally com- 
plicated pulsators, with closely-spaced periodicities that result in 
beating on time scales of hours to days or even longer. Furthermore, 
the intrinsic amplitude of certain pulsation fr equencies appea rs to 
change on the time scales of days to weeks ( Kilkenny 20 1(J). To 
make precise models of these stars, as well as to clearly determine 
rotation periods of days to months through multiplet splitting, we 
need to know which periodicities are intrinsic, and which show up 
in the Fourier transform because of amplitude (or phase) modula- 
tion. 

These effects are difficult to sort out with ground-based data 
because of diurnal aliasing, transparency variations and other atmo- 
spheric effects. As a result, disentangling the presence of closely- 
spaced (in frequency) periodicities from intrinsic amplitude varia- 
tion has been a challenge for many of the known sdBV stars. We 



are still uncertain whether intrinsic amplitude modulation is a ubiq- 
uitous feature of sdBV stars. Long, high duty cycle observations of 
sdBV stars should help us answer these questions. 

For their potential as probes of an important and poorly under- 
stood phase of stellar evolution, sdBV stars were included as tar- 
gets for the asteroseismology survey phase of the Kepler Mission. 
The Kepler scien ce goals, miss i on des ign, and overall perfor mance 
are reviewed bv lBorucki etalj d2010t) and lKoch et"alJd2010l) . As- 
teroseismology for Kepler is being conducted through the Kepler 
Asteroseismic Science Consortium (KASCfl These data were ob- 
tained as part of the survey mode, where short-cadence targets are 
observed for 30 days each during the first year of the mission 

We observed a number of candidate sdB stars during this 
phase, and have discovered several pulsators. The target selection 
process for sdB stars, and r esults of the fir s t 6 mon ths of the sur- 
vey phase, are described in l0stensen et ail djZOlOb) (hereafter Pa- 
per I). In this paper, we present the discovery of a new pulsating 
sdB star, Kepler Input Catalog (KIC) 010139564, at RA(2000) = 
19 h 24 m 58=15, Dec. (2000) = 47°7'53'.'6, that shows short-period 
pulsations and also shows evidence of being a hybrid (DW Lyn) 
star. Our results suggest that many of the periodicities are quite sta- 
ble in frequency, amplitude, and phase. Some of the periodicities 
are not as yet fully resolved in frequency, but we anticipate that fur- 
ther data will be able to determine whether their complexities are 
caused by intrinsic pulsations that are closely spaced in frequency, 
or are caused by amplitude and/or phase variations. 

The V316Hya stars, such as KIC 010139564, are in gen- 
eral much rarer than the longer-period V1093Her stars. Through 
the first 3/4 of the Kepler asteroseismology survey phase, 
KIC 010139564 is the only V361 Hya star found in the Kepler field. 
For first results from the Kepler s urvey phase for the long period 
sdBV stars see iReed et all (2010) (Paper III), with asteroseismic 
determinations of the p roperties of one of those stars presented in 
van Grootel et al.l fcoich (Paper IV). First results f or long period 
sdB p ulsators in binary systems are presented by iKawaler et al.l 
fcoich (Paper V). A bright, eclipsing sdB+dM system with a pri- 
mary that is pulsating with both p- and q- modes was also d iscov- 
ered in the Kepler commissioning data (0ste nsen et al .1201 Pel) . 



2 OBSERVATIONS 

KIC 010139564 has a Kepler magnitude (Kp) of 16.13 magni- 
tude. The Kp magnitude is a broadband magnitude reflecting the 
colour response of the Kepler photometer, with a central wave- 
length roughly equivalent to t he J?-band and a passband effec- 
tively spanning 437 to 897 nm dKoch et al]|2010h . This is a spec- 
troscopically confirmed sdB star (Paper I) with T c ff = 32, 500 ± 
200K and log g = 5.81 ± 0.04, and with a very small surface he- 
lium abundance: log(He/H) = -2.8. The data were obtained between 
BJD 2455002 and 2455033 (20 June - 20 July, 2009). During these 
~ 30.5 days of observations, a total of 45,210 science images were 
scheduled. However, spacecraft-related events occasionally result 
in loss and subsequent reacquisition of fine guidance, which re- 
sults in individual or short sequences of images being unusable. In 
addition, cosmic ray interference and software glitches can result 
in longer down-times (as the spacecraft enters "safe mode") for 
which no data are available. During this run, 3,686 images were ei- 
ther not obtained or were unusable. The photometric data obtained 

1 The Kepler Asteroseismic Investigation at the time of this work was man- 
aged by R. Gilliland, T. Brown, J. Christensen-Dalsgaard, and H. Kjeldsen. 
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is a result of the pulsational modulation. This is apparent when ex- 
amining the lower panel of Fig. [TJ which shows a periodic varia- 
tion of approximately 175 seconds, with a peak-to-peak variation 
of 30 mma. 

Note that at the time of this writing the pipeline for reduc- 
ing Kepler short cadence data is still being fine-tuned and tested. 
As a result, it is possible that some frequencies (and likely that the 
pulsation amplitudes) could be affected by the current Kepler re- 
duction pipeline with differing results in subsequent data releases 
as the process is improved. 
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Figure 1. Light curve of KIC 010139564. The horizontal axis is time (in 
days in the upper panel, seconds in the lower panel), and the vertical axis 
shows the differential flux variation in mma (parts per thousand). The upper 
panel shows the complete light curve, while the lower panel shows a shorter 
segment of the light curve from the first day of observation, illustrating the 
pulsations that are easily seen in the time series. 



for this target are short cadence data with a time between points of 
58.8 5 s and a duty cycle (within the 58.85 s window) of 92 per- 
cent dGillilandetalB oiO). Because of the gaps during the run, and 
other small deviations in the sampling rate, we determined an ef- 
fective Nyquist frequency of 8496.20 ^iHz by the symmetry in the 
power spectrum at frequencies surrounding this value. 

Raw and corrected photon counts are provided by the data 
pipeline. The "corrected" flux values included a preliminary esti- 
mate of the contamination from a nearby (assumed constant) star 
that contributes light to the 4 arc second square pixel in the focal 
plane. This correction subtracts a constant flux, which increases the 
resulting fractional amplitude of variations for the target star. For 
this data set, the estimated contamination fraction is 0.6 (that is 
40 per cent of the flux from the target and 60 per cent from back- 
ground or foreground stars). At this time we do not have data on 
the contaminating source, but we recognize that it could perhaps be 
variable, and contribute to peaks in the power spectrum. We expect 
that the contamination estimate (and the properties of the source) 
will be determined in the future, but for this work we report results 
for the raw fluxes. 

Kepler data on this star were provided with times in barycen- 
tric corrected Julian days, and flux. Variations of the normalized 
"raw" flux from the mean (of 1) are shown in the light curve in 
Fig. [U Note the gap during the data that corresponds to a saf- 
ing event. Amplitudes are given as parts-per-thousand (or milli- 
modulation amplitude, mma), with 10 mma corresponding to 1.0 
per cent. Since some contamination by background sources is in- 
cluded (see earlier discussion) the intrinsic pulsation amplitude is 
greater than this. Future processing (with an accurate estimate of 
contamination) should provide a more realistic absolute amplitude. 
The short period variations are clearly seen in the light curve - in 
fact, much of the observed width of the light curve in the top panel 



3 TIME-SERIES ANALYSIS 

The beauty of data from the Kepler spacecraft is that we get nearly 
continuous data acquired from a single instrument in a uniform 
manner. As such, analysis can proceed in a straightforward way. 
We begin with a standard Fourier transform / temporal spectrum of 
the entire light curve up to the Nyquist frequency of 8496.2 /jHz. 
The main frequency range showing oscillations is from 5000 to 
7000 /jHz. This region is shown in the top panel of Fig. [2] 

One point of concern for short period sdBV stars is that the 
sampling rate places the Nyquist frequency at 8496 [iHz, meaning 
that periods shorter than 117.7 seconds lie beyond the Nyquist fre- 
quency. This can be an issue a s a many sdBV stars h ave pulsation 
frequencies beyond that range ( 0stensen et al. 2010a). Shorter pe- 
riods will be mirrored across the Nyquist frequency into the lower 
frequency domain, creating spurious frequencies of intrinsic vari- 
ability to the star. In addition, since the oscillation frequencies are 
relatively close to the Nyquist frequency, and data are obtained by 
integration over the short-cadence interval of 58.8 seconds, the am- 
plitudes that appear in the temporal spectrum are reduced below 
the intrinsic amplitude. This reduction factor is approximately 25 
per cent for the shortest period frequency, and 13 per cent for the 
periodicities near 170-180 seconds. 

Also, though the data we examined were Kepler short cadence 
(SC) data sampled every 58.84 seconds, the spacecraft also sup- 
plies long-cadence (LC) data with intervals of 30 times the SC rate. 
Some part of the process by which the LC data are output affects 
the SC data by creating arte fact peaks at multiple s of the LC fre- 
quency of n x 566.427 pHz dGilliland et~aT1l2010l) . This artefact is 
not well characterized, but it is often most noticeable at higher har- 
monics that can impact the analysis of the short period sdBV stars. 
We examined data on many stars obtained during this same observ- 
ing period as KIC 010139564, and detected common artefact peaks 
at 3965, 4531, 5098, and 5664 /xHz. Therefore, we do not further 
consider periodicities at these frequencies in KIC 010139564. 



3.1 Fitting periodicities revealed by the temporal spectrum 

W e estimate the significance of peaks in the temporal spectrum as 
in iBregerl |l993): in frequency ranges near the range of interest, 
but that do not contain obvious significant peaks, we find the mean 
value of the amplitude spectrum (the square root of the power). In 
this case, in the frequency range of interest, the mean value of the 
temporal spectrum is approximately 37 parts per million, or an am- 
plitude of 0.037 mma. Defining that mean as a, we adopt an am- 
plitude threshold for significance of 4-cr. We note that this is just 
one option in determining a significance threshold; a more conser- 
vative approach is to define a as the mean amplitude of peaks in the 
temporal spectrum, since it is the peaks that are used in identifying 
candidate periodicities. However, here we will use the mean of the 
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temporal spectrum, rather than the peaks, to remain consistent with 
prior work. We therefore adopt a 4 a detection limit of 0.15 mma 
in the range of interest in KIC 010139564. 

Beginning with the highest peaks in the temporal spectrum 
we successively fit each periodicity (via non-linear least squares) 
in the time domain, and then removed each in the time domain, 
progressing down to the 4-a detection limit. In this procedure, we 
did not fit peaks that were within the frequency resolution limit of 
previously removed periodicities, nor did we fit peaks that were not 
visible in the original power spectrum. The frequency resolution of 
the data is formally 0.38pHz, which corresponds to the reciprocal 
of the time span of the data. 

The prewhitening procedure is illustrated in Fig. [2] In this fig- 
ure, the top panel shows the original temporal spectrum covering 
the region where most of the periodicities were found. The middle 
panels show the residuals after removal of seven and then nine high- 
amplitude frequencies, and the bottom panel has 22 frequencies re- 
moved (including two frequencies outside of the illustrated range, 
and one at the position of the instrumental artefact at 5098 /^Hz. 
The dashed (red) line indicates an artefact frequency mentioned 
above, the solid horizontal (blue) line indicates the 4-a detection 
limit, and the arrows indicate identified and fit periodicities. The 
arrows are color-coded (in the electronic version, shades of gray 
in the printed version) to indicate the stages of the fitting process. 
The magenta (lightest) arrows indicate the seven frequencies which 
were most readily resolved and fit, blue indicates other frequencies 
which were apparent after the prewhitening of the initial seven pe- 
riodicities. Some of these are extremely close to the initial seven, 
but others are quite well separated. 

In this star, there are at least four periodicities with ampli- 
tudes which would make them easily detectable from the ground. 
However, the bulk of the pulsations have amplitudes below usual 
single-site Earth-based detection limits for a 16th magnitude star 
observed with a l-2m class telescope. The solution to our fit, list- 
ing frequencies, periods, and amplitudes is provided in Table[T] 

While the pulsations span over 1400 fiHz, most frequencies 
are within a range of approximately 1000 /jHz. There are crowded 
regions where the periodicities are closely spaced or even unre- 
solved at the current resolution of 0.38 /jHz. The crowding is il- 
lustrated in Fig. [3] As in Fig. [2] the horizontal (blue) line is the 
4-a detection limit and the arrows indicate the positions of various 
periodicities that we included in the fit. Each vertical panel shows 
a 10 fjHz region with the top panels showing the original FT, the 
middle panels are prewhitened by the highest amplitude frequency 
in the region. The dashed (magenta) line indicates the original FT 
in this panel and illustrates the difficulties with removing additional 
frequencies. The bottom panels show the regions after prewhiten- 
ing at the indicated frequencies. The bottom panels are plotted on 
the same vertical scale, but the others have scales to optimize the 
visibility of peaks. Despite having 30.5 days of nearly continuous 
data, there remain at least three areas of unresolved periodicities. 
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Figure 2. Temporal spectrum of KIC 010139564. Top panel is the temporal 
spectrum of KIC 010139564 (the highest peak, at 10 mma, is truncated); the 
second panel is for the data with seven periodicities (indicated by magenta 
arrows) removed via prewhitening. Two additional periodicities (indicated 
in blue) were removed for the third panel. The bottom panel is the temporal 
spectrum of the residuals after all frequencies in Table 1 have been removed. 
The artefact at 5097 fiHz is indicated by a dashed (red) vertical line. The 
horizontal blue line shows the 4-ct detection limit. Note that the vertical 
scale changes with each panel. 



tion peaks which can have high frequencies 0. In this star, the first 
harmonic of f8 should be at a frequency of 11520.26 /jHz, which 
is 3024.26 /jMz above the Nyquist frequency. Thus an alias of the 
harmonic would appear at fa y — 3024.26 /iHz, or 5471.94 /jHz. 
There is indeed a peak in Tabl^T]at that frequency. Thus we identify 
f6 as the (aliased) first harmonic of f8. 

Similarly, high-amplitude sdB pulsators sometimes show pe- 
riodicities at linear combinations of strong pulsation modes. For 
the periodicities in this star, frequency sums will also be above the 
Nyquist frequency. We see one example of this (noted in TableQJ: 
the sum of f8 and f9 lies above the Nyquist, and its alias peak lies 
precisely (well within our errors) where f5 is found. Thus we con- 
clude that f5 and f6 result from non-linear behavior of f8 and f9. 
Similar ties to the dominant f8 mode link f3 and f 10 (i.e. the alias 
of f3 + f8 lies at precisely f 10) and f2 and fl 1 (the alias of f2 + f8 
lies within 0.11 pHz of fll). Because of the rather low Nyquist 



3.2 Combination frequencies, harmonics, and aliasing across 
the Nyquist frequency 

Many pulsating sdB stars show frequencies that are higher than 
the Nyquist frequency present in the Kepler SC data. Furthermore, 
large-amplitude pulsators often show harmonics of the main pulsa- 



2 For example the sdB pulsator PG 1325+101 ISilvotti et all |2006| ; 
ICharpinet et alj |2006|) has T e g and logg that is very similar to 
KIC 010139564 and is useful comparison. PG 1325+101 has a dominant 
pulsation mode at 7255.5 /iHz and its harmonic at 14,511 /iHz. It also 
shows other periodicities ranging from 5960 to 10,545/iHz. 
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Table 1. Periodicities in KIC 010139564 Formal least-squares errors in parentheses. 





Frequency [/iHz] 


Period [s] 


Amplitude [mm a] 


Qi-nl-v-il-itif fart 

oiaDiiiiy test 


Comment 


fO 


315.960 (30) 


3164.96 (30) 


0.435 (67) 






fi 


5286.046 (30) 


189.1773 (07) 


0.698 (66) 






t'2 


5287.298 (08) 


189.1325 (03) 


1.847(68) 


* 




f3 


5287.798 (12) 


189.1146 (04) 


1.173 (65) 




alias of f8+fl0 ? 


f4 


5413.385 (40) 


184.7273 (13) 


0.329 (68) 






f5 


5471.184 (76) 


182.7758 (25) 


0.187 (65) 




alias of f8+f9 


f6 


5471.930(46) 


182.7509 (15) 


0.297 (64) 




alias of 2xf8 


17 


5472.864 (03) 


182.7197 (01) 


3.891 (66) 


* 




f8 


5760.229 (01) 


173.60420 (04) 


9.998 (64) 


* 




f9 


5760.980 (04) 


173.5816 (01) 


3.314(64) 






f 10 


5944.369 (22) 


168.2264 (06) 


0.623 (66) 


* 




til 


5944.982 (39) 


168.2091 (11) 


0.349 (66) 




alias of f2+f8 ? 


112 


6001.343 (14) 


166.6294 (04) 


0.956 (66) 


* 




fl3 


6057.601 (55) 


165.0818 (15) 


0.236 (66) 






fl4 


6076.316(80) 


164.5734 (22) 


0.162 (64) 






fl5 


6172.424 (17) 


162.0109 (04) 


0.762 (66) 






116 


6234.673 (14) 


160.3933 (04) 


0.941 (67) 


* 




m 


6315.043 (77) 


158.3520 (19) 


0.169 (64) 






fl8 


6758.208 (90) 


147.9682 (20) 


0.145 (64) 






f 19 


7633.647 (75) 


130.9990 (13) 


0.175 (65) 







frequency for this star, and the complex structure of the temporal 
spectrum in the region of f2 and f3, it is difficult to determine, with 
confidence, which periodicities represent oscillation modes, which 
are combination frequencies, and which might be aliases of "true" 
peaks at higher frequencies. We hope that continued Kepler obser- 
vations will resolve some of these ambiguities, but higher-cadence 
ground-based photometry will also be needed. 

To summarize, we find 20 individual periodicities through a 
temporal spectrum / least squares fitting procedure, with 19 dis- 
tributed nonuniformly across the range of frequencies expected for 
an short-period sdB V star. Of those, at least two and perhaps four or 
more represent harmonics and linear combinations of the stronger 
periodicities, aliased across the Nyquist frequency. We find a single 
low-frequency periodicity at 316/iHz that is within the frequency 
range seen in the long-period sdBV stars. Thus, KIC 010139654 
appears to be a hybrid sdB pulsator (though there is a small chance 
that this periodicity could be caused by the contaminating source). 
In the next section, we examine the periodicities further to see if 
we can determine if they are all independent oscillation modes, or 
if some of the closely-spaced peaks are the signature of slow am- 
plitude and/or phase modulation. 



4 DISCUSSION 

4.1 Constraints on the pulsation modes 

For successful asteroseismic analysis, the observations need to pro- 
vide determination of (or useful constraints on) the modes of os- 
cillation represented by the observed periodicities. This includes a 
determination of the degree / of the oscillation mode. The mode 
frequencies themselves are described by spherical harmonics de- 
termined by the quantum numbers I, the mode order n, and (if ro- 
tation is present) the azimuthal quantum number m. When rotation 
is present, it can readily reveal the value of / if all 21 + 1 multiplets 
are present. From a stellar seismology perspective, this was demon- 
strated quite clearly in early gro und-based networ ked observations 
of pulsating white dwarfs (i.e. IWinget etal.lll99lh but has been 



used extensively for other c ompact pulsators as well (i.e. lReed et al.l 
l2004al ;lo' Toole et al .120041) . Time-series photometry rarely uncov- 
ers modes with £ > 2, so the multiplet structure should be well 
defined - unless the rotational splitting is larger than the spacing 
between modes with successive n, as can be the case in 8 Scuti 
stars (i.e. lTempleton et aO 19971) . 

In sdB stars, we expect that th e ir ro tation rates are slow 
enough jHeber, Reid. & Werned [l999l. 2000) that such overlaps 
should not be a p roblem jVan Grootel et al" I l2008t) except in 
some extreme cases jKawaler & Hostlerll2005ll For slow rotators, 
rotationally-split multiplets should be nearly equally spaced in fre- 
quency. Yet such multiplets are seldom seen in sdBV stars. This 
is where Kepler data should be most useful, as it should clearly 
resolve rotationally split multiplets in sdB stars even if they have 
rotation periods of months. However, data from the survey phase 
spanned only 30.5 days and so to resolve any multiplets, the rota- 
tion period would need to be shorter than that. 

We find no obvious multiplets in the periodicities of 
KIC 010139564. With a formal resolution of 0.38 /xHz, the small- 
est barely-resolved spacings that we find are 0.5 to 0.6 /^Hz. Hence 
we can state that the data are consistent with a rotation period cor- 
responding to P rot about 20 days or longer (1 / P m 0.6 /iHz). That 
said, the amplitude of m / multiplets does depend on the inclina- 
tion of the rotation axis to the line of sight, so if the axis is pointing 
towards the Earth the splitting would be at very small amplitude 
compared to the central peak. There are regions where the pulsa- 
tions are clearly unresolved and these regions are the most likely 
places to discern multiplets. We will need to wait for additional 
Kepler data to see if these regions resolve into multiplets. 

Without mode identifications via multiplets, we can, as a weak 
overall constraint, appeal to another property of sdB oscillations 
- the density of periodicities within a frequency range. Matching 
the observed T e g and log g with a stellar model allows us to es- 
timate the number of available mode frequencies within the ob- 
served range. In resolved sdBV stars, we sometimes observe many 
more pulsation frequencies than I = 0, 1, and 2 can provide, in- 
dependent of the number of inferred m 7^ modes. If higher 
I modes contribute to the observed number of periodicities, they 
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Figure 3. Enlarged sections of the temporal spectrum of KIC 010139564. The top panels are original temporal 
spectrum, and the middle and bottom panels are prewhitened by the indicated frequencies. The 4-<r detection 
limit is indicated by the horizontal (blue) line and the dashed (blue) line in the middle panels is the original FT, 
to show how the high-amplitude peaks affect the others in that region. Note that the vertical scale changes with 
each panel. 



must have a higher intrinsic amplitude than the low er I modes be- 
cause of the large degree of geometric cance llation jDziembowskil 
1 1971 : ICharpinet et al]2005l : lReed et al.ll2005[) . 

Models appropri ate for KIC 010139654, such as those in 
ICharpinet etafl 12006) forPG 1325+101 (which has similar values 
for Tog and log <?), show roughly one overtone (n) per I degree per 
850 /iHz. Under normal observing conditions, this would supply 
three to four modes per 1000 pHz (I = 0, 1, 2) without invoking 
higher t values and if no multiplet structure is present. 

For KIC 010139564, the peaks fl to fl7 in Table [j] span 
slightly more than 1000 /jHz. If the closely-spaced periodicities are 
considered as multiplets of one overtone, then 10 modes remain, 
which cannot be accommodated by low I modes 

The amplitudes of the Kepler periodicities in Table Q] span 
nearly a factor of 70, so that it is possible that the lower- amplitude 
periodicities might have £ > 2. If further Kepler observations span 



3 If we allow for 1 = 3 and £ = 4 modes, recognizing that I = 4 
modes suffer less geometric cancelation than I = 3 modes for most m 
values ICharpinet et al. 2005; Reed et al. 2005), the number of modes per 
1000 /^Hz rises to 6. 



one year, and the modes remain stable, the noise level will be re- 
duced by a factor of 3 or more, and if higher I modes are present 
we might hope to see them. 

4.2 Amplitude and phase stability 

If pulsating sdB stars are observed over an extended time period, it 
is common to detect amplit ude variability in many, if not all, of the 
pulsation frequen cies (e.g. lO'Toole.et al.ll2002l ; lReed et al.f 2007b; 
IZhou et aljibOrj) . Variations on time scales of days to weeks some- 
times appear in clearly re solved pulsation sp ectra where mode beat- 
ing cannot be the cause dReed et alf2006ah . Given the limitations 
of ground-based observations, in most cases long-term variability 
cannot be securely attributed to beating between pulsations, and 
could also result from intrinsic amplitude or phase modulations. 

As previously mentioned, a true strength of Kepler data is that 
we have nearly continuous, single-instrument data which can be ex- 
amined for amplitude and phase stability. We broke the data up into 
six 5-day spans, and computed an averaged temporal spectrum of 
these six shorter runs to obtain estimates of the position of the cen- 
troid of pulsation "power" at each of 7 cleanly separated regions in 
the averaged temporal spectrum. We then used those 7 frequencies 
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Figure 4. Pulsation phases for 7 well-separated strong pulsations seen in 
KIC 010139564. Here, we used a least-squares fit to the dominant frequency 
in each band (with fixed frequency) and fit 2-day subsets of the data to 
obtain the pulsation phases and amplitudes (see Fig. [5). The phase is ex- 
pressed in terms of the time of the first maximum (in seconds) after the start 
of the observations. The three highest-frequency periodicities (bottom pan- 
els) show no phase variation over the course of the run, but the four lower 
frequency periodicities show apparent variations. 



(flagged with asterisks in Table [T} as fixed quantities to examine 
them for phase and amplitude stability over two-day spans of data. 
We used a fixed value for each of these 7 frequencies (see Table[T} 
and then did a non-linear least squares (NLLS) fit to the data in 2- 
day intervals to find the best-fitting values for their amplitudes and 
phases. The results are shown in Fig.|4]and Fig. [5] 

Breaking the data into two-day chunks reduces the frequency 
resolution, but it allows us to track the amplitude and phase of the 
variation with time. If the amplitude is modulated because of beat- 
ing between two coherent modes, then the phase of the light curve 
should flip sign at the time that the amplitude reaches a minimum 
(or maximum). Both phase and amplitude curves should look sinu- 
soidal as well, as if there are only two modes beating. If the ampli- 
tude or phase curves show non-sinusoidal variations on time scales 
shorter than the run length, and the extrema in the amplitude vari- 
ations do not occur at zero-crossings of the phase variations, then 
simple beating between two or three modes is not taking place. 

From Fig. [4] we see that the three higher frequencies examined 
clearly demonstrate phase stability over these time-scales. Their 
amplitudes are also stable (Fig. [5}. Over the course of Kepler's 
three and a half year mission, these frequencies will be useful for 
searching for secular period changes caused by the evolution of the 



Figure 5. Pulsation amplitudes (in mma) for the same 7 well-separated 
strong pulsations used in Fig. [4] computed using the same procedure. The 
three highest-frequency periodicities (bottom panels) show no obvious am- 
plitude variation over the course of the run, but the four lower frequency 
periodicities show apparent variations. 



stars, and ma y be useful for prec ise timing to search for sub-stellar 
companions dSilvotti et al.ll2007h . 

The four lower frequency periodicities (top 4 panels) show 
phase and amplitude variations over the 30.5 days of data. Of those, 
two (5472.9 /iHz, f7 and 5760.2 fiHz, f8) vary on time scales 
shorter than the length of the observing run. Their phase and ampli- 
tude variations are out of phase by about 90 degrees, as one would 
expect from beating between closely spaced stable oscillation fre- 
quencies. Thus we conclude that these two frequencies appear to be 
(nearly) resolved into two or three stable periodicities. The largest 
amplitude periodicity (5760.2 /iHz, f8) shows sinusoidal amplitude 
and phase variations, and is cleanly resolved in the temporal spec- 
trum analysis (see Table Q] and Fig.|2j with separations of approx- 
imately 0.75 /jHz. The 5472.9 /iHz region can also be fit by three 
periodicities spanning 1 .7 /iHz, but not with the same level of con- 
fidence - note also that two of those may be linear combination 
frequencies. 

The peak at 5944.4 /jHz (flO) also shows variation in ampli- 
tude in Fig. [5] but the phase variation shows a less regular pattern 
that does not match the expected phasing. In part this may be a re- 
sult of the pulsation amplitude reaching small values (i.e. around 
day 20). The power spectrum analysis suggests that this is a pair of 
peaks separated by approximately 0.61 /iHz, but it is not well re- 
solved, and again may include a periodicity at a linear combination 
frequency. 

The spacings hinted at above, based not on the least-squares 
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fits in Table [T] but on the beating pattern, are all around 0.6 to 
0.8 i-iHz - roughly twice the frequency resolution of the data. If 
we attribute these splittings to rotation, then we can suggest a ten- 
tative value for the rotation period of approximately 16-20 days for 
KIC 010139564. Clearly a longer run is needed to firmly resolve 
these peaks and confirm that they are indeed multiplets. 

Given the current data, the best case we find for a periodic- 
ity that has an intrinsic variation is the region around 5287.3 /jHz. 
This shows phase variation on a time scale longer than the span of 
observation. The amplitude varies quasiperiodically with a period 
roughly 0.75 times the observation span, but on a different (shorter) 
time scale than the phase variation. In the temporal spectrum, this 
is a region where NLLS fitting failed, and residual power remained 
after prewhitening. Given the phase and amplitude variations seen 
in the data, we suspect that this mode shows intrinsic amplitude and 
phase instability, unlike the other strong periodicities in this star. 

4.3 The low-frequency oscillation 

The periodicity labelled fO in Table Q] is notable for being in the 
frequency range observed in the V1093 Her stars, and near the low- 
frequency periodicities in the hybrid (DW Lyn) pulsators. This sin- 
gle periodicity, with an amplitude of 0.45 mma is the smallest- 
amplitude variation seen in this frequency range for hybrid subd- 
warf B pulsators. In the other members of this class, the amplitude 
ratio between the short and long period modes is in the range from 
5 to 20 or more. In this star, the ratio is a factor of 23, which is 
consistent with that seen in the known hybrid pulsators. 

Fig.[6]shows the temporal spectrum in the low frequency range 
surrounding fO, along with the time dependence of the pulsation 
amplitude and phase (computed in the same manner as for the short 
period modes). A few things are notable about this mode. First, it 
seems to be a single oscillation mode - fitting a single sinusoid re- 
moves the temporal spectrum peak and leaves no residuals that can 
be distinguished from noise. It remains visible through the entire 
30.5 day interval, and maintains phase (within the errors). How- 
ever, there is a possible trend towards decreasing amplitude over 
most of the observing interval. With further monitoring, we hope 
to determine if this mode will remain visible. 

Interestingly, with the spectroscopically determined T e g of 
32,500K, KIC 010139564 is significantly hotter than the other hy- 
brid pulsators. Finally, we add that there is a slim chance that this 
periodicity may be contributed by the contaminating source - if 
Nature is being perverse. 



5 CONCLUSIONS AND FUTURE WORK 

We present results from a study of the first Kepler data on a newly- 
discovered pulsating subdwarf B star that shows mostly short- 
period pulsations characteristic of the V361 Hya stars, along with 
a single oscillation at a frequency characterstic of the VI 093 Her 
stars. Thus KIC 010139564 may be a member of the hybrid sdB 
pulsator class. The Kepler data have realized the potential to ob- 
tain a nearly uninterrupted time series of photometric data for this 
type of pulsator that spans over 15,000 pulsation cycles. Previous 
ground-based efforts have demonstrated that such data are essential 
for overcoming aliasing issues, caused by gaps and the necessity of 
using multiple instruments on varied apertured telescopes from the 
ground. The science questions raised by ground-based attempts at 
photometry of these stars can be addressed directly with Kepler 
data. 



These data have allowed us to confidently detect 21 period- 
icities, a few of which may be linear combination frequencies. At 
least one of the main periodicities remains unresolved, but it is not 
yet clear whether we are seeing beating or stochasitic variations. 
Residual power indicates the presence of possible pulsations below 
the current detection limit. As such, with longer duration Kepler 
observations, we can anticipate finding more periodicities at even 
lower amplitudes. Kepler can place much more stringent upper lim- 
its on pulsation amplitudes in these stars than we can achieve from 
ground-based data, so we can also hope to determine if ground- 
based observations of these stars that reveal simple pulsators have 
been limited by a higher detection threshold. We also anticipate 
resolving rotationally-split multiplets, which will add stringent ob- 
servational constraints to theoretical models. 

The data presented here confirms the potential of Kepler for 
sdB asteroseismology. KIC 010139564 has a frequency spectrum 
that is as rich as (or richer than) possible for low-degree modes 
only, and there is still more to come. We anticipate further, longer- 
duration observation during the second year of operations, when 
Kepler will no longer be in survey mode. 
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